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Abstract 

The shapes of cell membranes are largely regulated by membrane associated, curvature active, proteins. We use 
a numerical model of the membrane with elongated membrane inclusions, recently developed by us, which posses 
spontaneous directional curvatures that could be different along and perpendicular to its long axis. We show 
that, due to membrane mediated interactions these curvature inducing membrane nematogens can oligomerize 
spontaneously, even at low concentrations, and change the local shape of the membrane. We demonstrate that for 
a large group of such inclusions, where the two spontaneous curvatures have equal sign, the tubular conformation 
and sometime the sheet conformation of the membrane are the common equilibrium shapes. We elucidate the 
factors necessary for the formation of these protem lattices. Furthermore, the elastic properties of the tubes, like 
their compressional stiffness and persistence length are calculated. Finally, we discuss the possible role of nematic 
disclination in capping and branching of the tubular membranes. 
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I. INTRODUCTION 



Membrane shape deformations are key phenomena in a multitude of cellular processes, including pro- 
tein sorting, protein transport, organelle biogenesis and signaling. In the last decade a profusion of 
regulatory proteins facilitating such shape changes of the cellular membranes has been unravelled, with 
the BAR protein superfamily [1], the Pexll family [2] and coat proteins [3] as notable examples. The 
possibility of such mechanisms has long been anticipated in the biophysical literature[4, 5]. However 
the experimental and theoretical difficulties involved have hampered the establishment of a quantitative 
basis for interpreting such phenomena in cell biology. Recently, we had overcome one such obstacle by 
the establishment of a computer simulation technique to study how the cooperative effects of membrane 
inclusions, imposing a curvature along the direction of its orientation, remodels vesicular membranes[6]. 

In this work we aim at describing, from a theoretical point of view, the effect of a large group of these 
membrane curving proteins, which can be considered as effectively elongated objects in the plane of the 
membrane. We consider inclusions with approximate 7r-symmetry, i.e. the protein can be considered 
as essentially indistinguishable from its form rotated by 180° around the protein center in the plane 
of the membrane. The membrane inclusions we consider, has thus some similarity with nematogens 
in 3-D nematic liquid crystals. However, they are embedded in a membrane and may couple to its 
geometry, and it is only the part of the protein in contact with the membrane, which will be subject 
to these symmetry requirement. Therefore, we cannot consider these membrane inclusions as simple 
liquid crystal nematogens restricted to Euclidean two dimensional surfaces. In the following we will refer 
to such membrane inclusions as membrane nematogens. Large groups of membrane curving proteins 
fall into this category of membrane nematogens. An example is the BAR proteins(proteins containing 
both BAR domains and/or N-terminal helices), where both the N-terminal amphipatic helices and the 
banana-shaped, positively charged, dimeric interface with the membrane, induces directional curvature 
[7-12]. The caveolin protein family [13], which form dimers and are bound to the membrane by a pair 
of hairpins and the reticulon, DPI and Yoplp involved in the formation of smooth ER [15, 16], and 
are anchored to the membrane by two similar hairpins are also examples. The cell biology literature 
has provided good evidence for that the insertion of amphipathic helical peptide sequences not only 
provide a binding mechanism, but also gives rise to local modulation of the membrane curvature [17, 18]. 
More solid, quantitative support for this conjecture is given from biophysical experiments [19] and theory 
[20] . Furthermore, biophysical studies has demonstrated that curvature active membrane inclusions 
have dramatic effects on the cooperative behavior with a complex interplay between lateral ordering and 
membrane shape. However, the detailed mechanisms leading to the specific complex membrane-protein 
structures have not been characterized. This work will elucidate some aspects of these mechanisms for 
the membrane nematogens. 

Some of the key processes involved in the structural organization of membrane nematogens described in 
the cell biology literature, can be categorized as follows: (1) the aggregation of the nematogens - the pro- 
cess where membrane proteins upon activation and/or binding to the membrane spontaneously aggregate 
and form functional cluster of proteins in the membrane [9-11, 14], (2) Tubulation of membranes, where 
the aggregate and the membrane develop tube-like membrane structures (e.g. sorting endosomes[21, 22] 
and Mitochondrial outer membrane [23], formation of T-tubules in Drosophila[24:]) and (3) The formation 
of protein lattices, often characterized by helical arrangement of the proteins spiraling around the tubular 
membrane, e.g. for dynamin [25, 26] or caveolin [27]. 

In this work, we will demonstrate by numerical analysis of a possible physical model, which captures 
the membrane conformations and the organization of in-plane nematogens, that the above mentioned 
processes directly results from the cooperative thermodynamic behavior of the nematogens coupled to 
the flexible membrane. Also, we will discuss aspects of the stability of membrane tubes and the formation 
of the edges for membrane sheets. Our model gives a coarse description of the membrane, which capture 
properties of the membrane which are essential for its large scale organization. Despite the simplicity of 
the model, the parameter space is too large for a comprehensive discussion of it's phase behavior. Rather, 
we will focus on some generic features of the model which may well give a framework for interpreting the 
experimental observations of cellular membrane morphogenesis. Previously, protein induced membrane 
tube formation has been considered by a phenomenological model involving scalar fields [28], and the 
coupling between membranes and inclusions with directional curvature was modelled in [29-32] . 

The paper is organized as follows: In Section II the physical model of the interacting system of mem- 
brane and membrane nematogens are presented, while details about the numerical analysis is given in 
Supplementary Materials. Section III on Results and Discusion present some generic properties of 
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the model and discuss their possible relevance to experimental results. In section IIIA-IIIC the aggrega- 
tion of proteins and membrane domain formation, membrane tabulation and formation of protein lattices 
are described in the framework of the model. Section III D discusses the elastic properties of membrane 
tubes and their relevance to observable effects. Much of the characterization of the elasticity of protein 
lattices is based on a continuum version of the model discussed Supplementary Materials. Section 
III E discusses mechanisms of closing, capping and branching of membrane tubules and the possible role 
of nematic point defects. Section IIIF describes aspects of the stability of membrane tubules with more 
membrane curvature components. In section III G the interplay between sheet and tubule formation is 
described and possible implications for cell organel morphology is discussed. Some perspective on the 
modeling of membrane morphogenesis is given in Conclusion, Section IV. We will in this work specialize 
to properties of membranes with inclusions which posses directional curvatures of equal sign. We will 
consider cases with different signs of the directional curvatures in a seperate publication. 

II. MODEL 

The modeling of the effects of in-plane nematogens on membrane structure, will in this work, be 
treated with a discretized description of the surface as a randomly triangulated mesh. A continuous 
surface conformation is approximated by a collection of triangles glued together to form a closed surface 
of well defined topology. A triangulated surface, with spherical topology, thus consist of N vertices 
connected by = 3(7V — 2) links, which enclose Nt — 2{N — 2) triangles. Each vertex v is assigned 
a position Xy. This tesselatations of the surface form the basis for a coarse grained description of the 
membrane, where only the gross features of the structure and interactions are important. 




(a) (b) (c) 

FIG. 1: (a) A one ring triangulated patch around a vertex v. The shaded region represents the tangent plane 
at V and N{v) its corresponding normal. ci„ and C2u are the maximum and minimum principal curvatures, 
respectively, along principal directions ti^ and t2v (b) Illustration of the nematic field vector n defined on the 
tangent plane of vertex v. (c) A vesicle of spherical topology with spatially random surface nematogens. 

The triangulation and the vertex position form together a discretized surface, a patch of which is given 
in fig:??. The geometry of the continuous surface, which is approximated by the triangulated surface, 
can now be characterized by a number of surface quantifiers, e.g. the curvature tensor, the principal 
directions {tiv&^t2v), the corresponding principal curvatures (cii, & C2v) and surface normal, N{v), at 
each vertex v. The details can be found in [6]. The discretized Helfrich's free energy [34] can then be 
evaluated as 

v—l 

where Hy = (ci„ -|-C2«)/2 the mean curvature at vertex v and Ay is the area of the surface patch occupied 
by the triangles adjacent to vertex v. k is the bending rigidity of the membrane. Furthermore we are in 
a position to calculate the directional curvatures along and perpendicular to a unit vector h along the 
surface by use of Gauss formula: 

Hy^W = ClyC0S^ipy+C2ySin'^^Py, 
Hy^l_ = Cly Sin^ Ify + C2yCOs'^ ify, (2) 
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where is the angle between fi and the principal direction ti^. Such an orientational spontaneous 
curvature may be induced by a membrane nematogen with an orientation in the plane of the membrane 
given by h. In addition to the interaction with the membrane, nematogens may tend to orient along 
each other at close proximity due to steric, electrostatic and dispersion interactions [35]. In the present 
study, we focus only on the two dimensional orientational interactions promoted by the underlying, non- 
planar, fluctuating membrane [36-40]. The 7r-symmetry of the individual nematogens dictates that the 
simplest form of their self interaction should be of the type cos^(0^u) and sin^(0t,„), where 6^ is the angle 
between fi^ and at neighboring vertices. We choose to represent the interactions between membrane 
nematogens by an extension of the well-established Lebwohl-Lasher model of nematic ordering in presence 
of vacancies, here implemented on a triangulated surface model of a membrane. 

The nearest neighbor interaction between the nematogens is composed of an isotropic component 
represented by an interaction strength J and an anistropic (quadrupolar) correction measured by the 
interaction constant eLL. The total interaction between the membrane nematogens thus takes the form 

Hfield = Il{-f-eLL(^ COS2(0,„) - Ml (3) 



where the sum is over nearest neighbour vertices. = 0, 1 is an occupation variable, which is unity if 
vertex "w" is occupied by a nematogen and zero if otherwise. The calculation of the d^u is non-trivial, 
since the angle between spatially separated nematogens are measured after the parallel transport of vec- 
tors along the curved surface [6]. With this measure of the angular differences, Eq.(3) models the in-plane 
interactions of the nematogens mediated by the membrane. The direct distance dependent interactions 
through the cytosol is not taken into account in this model of membrane-protein conformations. Suf- 
ficiently large, positive ell favors in-plane ordering of the nematogens. The effect on the anisotropic 
elasticity of the membrane due to the nematogens, in analogy with the discretized Helfrich free energy, 
takes the form [6]: 



and are the spontaneous curvatures along h and n , while Ky and are the corresponding 
directional membrane bending elastic constants. 

Self avoidance of the discretized surface is ensured by imposing constraints on the neighboring vertex 
distance and on the dihedral angles between neighboring faces [6]. The equilibrium properties of the 
discretized surface can now be evaluated by standard Monte Carlo sampling of Boltzmann's probability 

distribution ^ exp ^— ^^[Hc + "Hfieid + ^nc]^ at fixed concentrations of the membrane nematogens. 
A general description of the implcmetation of such numerical models and further details about the 
simulations are given in [6]. 

Finally, we will make some considerations about length scales. The lattice model is a highly coarse- 
grained representation of the membrane, designed to capture the large length-scale properties of mem- 
branes with inclusions. Therefore, the triangulated surface represent a collection of membrane patches 
with a characteristic length scale. A natural choice of length scale is to identify a tether length with the 
lateral extension of a membrane inclusion. Some examples here are CIP4 F-BAR with a length of 22 
nm[41] or dynamin which extent about 25 nm [25]. 

The computer simulations of the discrete model provide us with insight into the nature of equilibrium 
configurations for a choice of model parameters. To complement the numerical simulations, it is useful 
to consider the corresponding continuum model in the limit of membrane nematogen with 100% surface 
coverage. It is an extension of Helfrich's bending free energy functional [42, 43] 



<^ dA|^Tr(Vn : Vn) + ^(2i?)2 



where Ka — ^-^^ll- In Supplementary Material is presented such an analysis of the mechanical 
properties of a tubular membrane with a protein coat and analytical expressions reflecting, tube radius, 
persistence length and protein organization are also derived. 
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III. RESULTS AND DISCUSSION 



In this section we will present some key aspects resulting from the coupling of membrane nematogen 
proteins to lipid membranes. It will both contain results from computer simulations of the aforementioned 
model, which are non-perturbative, along with theoretical analysis of the continuum model, of a more 
perturbative character to qualify the numerical finding. Throughout the discussion the parameter e^^ 
has a relatively high value (several ksT in a range where nematic ordering is favored). Furthermore, the 
implications of our results on the experimental systems in vivo and in vitro will be discussed. 

A. Aggregation and membrane domain formation of membrane nematogens 

A common feature of membrane nematogens is their strong tendency to self-associate, driven by the 
flexible geometry of the membrane - in this manuscript, we call this self associated structure to be an 
aggregate or a domain. Self association has been observed for a wide range of model parameters 
K±, and H^. All results presented in the following corresponds to system size with N ~ 2030 
vertices. When the fraction of nematogens (j)A — 0.3, eLL — 3 and J = (in units of ksT), the flexible 
membrane with curvature coupled to the nematic orientation, gives rise to co-existence of nematically 
ordered domains and the isotropic dilute phase, this is shown in figure-2(b). This is to be compared with 
the planar Lebwohl-Lasher model on a random triangular lattice, at the same concentration, where the 
isotropic phase is stable (see Supplementary Materials). Additional direct repulsive interactions J < —0.5 
between the membrane nematogens can reestablish the isotropc phase, which is shown in Fig. 2(a). 
The aggregation of membrane nematogens cause shape deformation of the whole membrane with the 
collective involvement of all the degrees of freedoms, lateral orientation, lateral position and membrane 
conformation. In general the lateral domain formation depends on all the involved parameters - e.g. 
increasing J promotes the aggregation and can change the aggregate shape as shown in Fig. (2)(c). 




-l<J<-0.5 -0.5<J<0.5 0.5<J<1 



FIG. 2: Equilibrium membrane conformations, with (I>a ~ 0.3, k = 20, = 5 and H^^ = 0.5 and eLL = 3, for 
different range of J. 

The effect of concentration is shown in Fig. (3) for surface coverage in the range 4>a — 0.1 — 0.7, which 
display a series of complex shape deformations connected to different aggregate structures. More details 
will be discussed in section III G. 

The aggregation of membrane nematogens also has a temporal aspect. In Fig. (4) we have shown a 
Monte Carlo time series, for a membrane coverage of 10% nematogens, to illustrate the sequence of 
domain formation and membrane curvature induced changes leading to the equilibrium structure. The 
membrane nematogens in an initial randomly dispersed orientation assemble into smaller orientationally 
ordered domains mediating the final equilibrium structure. These ordered domains often appear as 
metastable configurations, which either disperse again due to lateral fiuctuations or they will eventually 
funnel into a equilibrium domain configuration. The Monte Carlo dynamics does not reflect the physical 
kinetics very well, but is useful in identifying kinetic paths connecting various metastable states that lead 
to the global minima [44, 45] . 

The aggregation of membrane inclusions mediated by membrane curvature deformations and fluctuations 
is not speciflc for nematogens, but is a more general phenomena for membrane curvature active com- 
ponents. It has been well understood in the framework of models for curvature instabilities [5, 19, 46], 
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(a) 0^ = 0.1 (b)(/.A = 0.2 (c)0A = O.3 (d)0A = O.4 




(e) (Pa = 0.5 (f) = 0.6 (g) cf^A = 0.7 



FIG. 3: Equilibrium configurations for varying composition, k — 10, = 5, k_l = 0, H^^ = 0.5, — 0, J — 0, 
4>A = 0.1 — 0.7 and tLL = 3. 




(a) Tmcs = (b) Tmcs = 1-6 x 10^ (c) Tmcs = 17.5 x 10^ 

FIG. 4: Aggregation of membrane inclusions for k = 10, Ky — 5, k± — 0, //y* — 0.5, — 0, J — 0, cj>A ~ 0.1 
and ELi = 3. Monte Carlo time series showing a) random initial configuration of membrane nematogens, b) 
intermediate state with multiple nematic domains and c) equilibrium conformation where all the small domains 
coarsen into a single patch. 

and has also been demonstrated that simple amphiphathic inclusions, e.g. antimicrobial peptides like 
Magainin or Melittin [19, 47] and viral membrane active proteins like NSB4 of Hcpctites C[48]. 

The self-association of these membrane components thus needs not to be facilitated by strong direct 
attractive interaction amongst them. The coupling to the membrane geometry provides additional indirect 
membrane conformation mediated attractive forces making them to slip into bound structures involving 
both the proteins and membrane curvature. However, the structure of the aggregates are dependent on 
the details of the molecular structure and the direct interactions. A general feature of these aggregates is 
that they appear as nematically ordered domains, where the nematogens form elongated oriented patches 
with well defined curvature characteristics, e.g. ridges or cylindrical rims. In the following we will in 
particular focus on the tube-like structures. 
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B. Tube formation 



The most prevalent equilibrium domain structure is the nematic tube, where the membrane protrude 
into a cylinder like configuration with the membrane nematogens forming a coat around the cylinder. 
Also for tube formation the overall interaction strength (J) between the membrane nematogens plays a 
secondary role. It's most pronounced effect is to widen the concentration range for tubulation and to 
enhance the line tension at the domain boundary, which can induce fission of tubes by narrowing the tube 
at the boundary of the domain, as in Fig(2(c)). The effect of concentration of membrane nematogens on 
the membrane tubulation phenomena is shown in Fig. (3(g)). For large concentrations of nematogens or 
increasing values of J the tubes are the characteristic equilibrium structures shown in Figs. (2,3). 

The radius of the equilibrium membrane tubes appears to be relatively well-defined. The radius of the 
tube with nematic order can be calculated on basis of the continuum model Eq.(5) for the chosen model 
parameters (see Supplementary Material): 



So, the radius f is set by the curvature elastic model parameters, involving the absolute value of the 
directional spontanous curvatures, modulated by the curvature elastic constants. It follows from Eq.(5) 
that the actual tube radius is somewhat larger than the inverse directional spontaneous curvatures and 
dependent on the relative strength of the elastic constants. 

In experimental systems the membrane tube dimensions can vary considerably with different types 
of proteins in the cell[l]. Membrane tubes formed in vitro by curvature active proteins also display a 
considerable variability in size. Frost et al. [41] have studied the effect of a number of mutants of CIP4 
F-BAR on liposomes. By mutations they find a big variations of tube diameters in the range of 50 to 



Membrane tubes induced by membrane inclusions are common phenomena in biological cells, both 
as more static structures like T-tubules of the muscle cells [24] or more temporal structures like sorting 
endosomes [49]. The examples shown in Fig. (2,3) corresponds to the cases where spontaneous curvatures 
are positive, like that induced when F- BAR-domain proteins bind to organelle membranes. However, if 
the proteins induce negative spontanous curvatures, as in I-BAR domain proteins, it gives rise to tubular 
invaginations as shown in Fig. (5) . 

As can be seen from Fig. (5), for proteins with large negative spontaneous directional curvatures, at low 
concentration {(j>A = 0.1 — 0.3), we obtain tubes growing into the interior of the vesicle. As 4>a increases, 
tubes disappear and saddle like regions appear. The inner tubes and saddle like regions coexist again for 
large concentrations (pA > 0.8. 

C. Protein lattices 

Membrane nematogens organize as nematically ordered domains and coat around the membrane to 
form tubes. Nematogens orient perpendicular to the tube axis when k± — 0, K|| ^ and 77" > 0. 

Similarly, Ky = 0, ^ and Hj_ ^ leads to an arrangement of the nematogens along the tube 
direction. For the common membrane nematogen both these parameters are non-vanishing. Such a case 
is shown in Fig. (6). 

The helical arrangement of the membrane nematogens at the tube surface can be easily understood 
considering that in general such arrangement will give rise to a global nematic ordering of the membrane 
nematogens (generalized spirals are the only geodesic curves on long tubes) and the radius is set by the 
elastic terms. The coupled expressions for the mean values for tube radius f and the angle ip between the 
tube direction and the nematogen orientation, for different regimes of the dimension less parameter tp: 




(5) 



100 nm. 
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(a) (l)A = 0.1 (b) = 0.2 (c) (Pa = 0.3 (d) 0a = 0.4 




(e) (pA = 0.5 (f) (Pa = 0.6 (g) (Pa = 0.7 



FIG. 5: Equilibrium configurations for vesicle with negative spontanous curvatures, k = 10, = 5, k± = 0, 
H?^ = -0.5, Hi = 0.0, J = and (f>A = 0.1 - 0.7. 



\2 



for V < 0, 



r = < 



The parameter ip is given by the model parameters as: 



K / 1 1 



(6) 



(7) 



similarly for the angle ip: 



r 



C0S^((/3) = < 



K± + KU 



for ^ 0, 



< V < 1 



(8) 



1 



-0 > 1 



So, in general both ip and r are set by the model parameters. A derivation of the above expressions are 
given in Supplementary Materials . 
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FIG. 6: Protein lattices. Modes of a tubular membrane at different state points, for k = 10 and eLL = 3. (a) 
Tubular conformation with {ip) — for k\\ = 5, k,± = 0, H'^^ = 0.4, (b) Spiral modes of the tube with (ip) — 
seen for Ky = 5, k± — 0, -ffij* — 0.4, and (c) rearrangement of nematics into spiral modes ((<p) ^ 0) when Ky = 5, 
K_L = 5, //? = 0.4, Hi = 0.25 



The spiral organization of the membrane coating proteins has now been observed for many tubular 
membrabe systems in vivo and in vitro, e.g for the F-BAR proteins[26, 41]. EM-tomographs of tubules of 
CIP4 F-BAR on liposomes [41] show a fairly dense packing arrangement in the helical tube. The average 
tube diameter is around 68 nm and the helical angle is about ip = 40°. Such arrangements are termed 
protein lattices in the cell biology literature. It is found that the helical angle f of the protein lattice with 
respect to the tube direction adjusts to the tube diameter such that the directional curvature is about the 
same. For a similar type of experiment with dynamin[50] the membrane tubes of radius r ~ 23nm with 
densely packed helical dynamin coat was observed with a helical pitch of about 15nm (corresponding to 
a helical angle (p — 80°). 

Our simulation results suggests that the spiral organization of the protein coat on the tube need not 
be a result of polymerization as often referred in the literature, but can be a self-assembly process of 
the curvature active proteins mediated by the membrane. Furthermore, the modelling suggest that these 
protein lattices are not conventional two dimensional lattice structures like polymerized membranes or 
graphene, but rather two dimensional nematic liquid crystalline structures. In the model there is no terms 
which can distinguish between a right or left turning helix, i.e. the helical arrangement is the result of 
a spontaneous symmetry breaking. However, the smallest chiral symmetry breaking contribution to the 
free energy can favor one of the helical orientations without having an effect on any other parameters. 

D. Thermal stability of membrane tubes 

While our model parameter determine the mean physical properties of the tubes , e.g. the radius, we 
expect the tubular membranes to display an elastic response to deformations in its shape and organization 
of the membrane nematogens. This can, for e.g., be reflected in the variation of the shape characteristics 
due to thermal fluctuations. For analysis of such deformations the continuum description of coated 
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membrane tubes are suitable and the details can be found in Supplementary Material. It is shown 
that in general the deviations in the orientation of the membrane nematogen and the tube radius are 
strongly correlated. The thermally induced fluctuations in the radius is found to be 

«^ = M: ^Lt^forO<^<l, (9) 

f2 47r K||K_L + (k|| + V- > V ; 

where f and Cp are respectively the equilibrium tube radius and nematic orientations and i\} = cos^ ip. 
We note that the relative variance in r has an upper limit "If^. With a typical range k ~ 20 — SOkeT 
this ratio in Eq.(9) is of the order 0.01. For CIP4 F-BAR, reconstituted on liposomes, cryo-tomography 
measurements give f — 33nm and {{Sr)'^)/f'^ ~ 0.01 [41]. If this observed variation in tube thickness 
is interpreted as frozen in thermal variations, it is in agreement with the the above theory. For rigid 
membranes with large k and/or large kj_ values we can consider the thermaly excited variations in r 
as small. Similarly, we can estimate the thermal fluctuations around cos^((^) for such a segment as, 

To our knowledge no experimental reports on the random variations in the helical angle has been given. 
A third type of deformation to consider is the bending of the tubes. It is shown in Supplementary 
Materials that when r is a constant along the tube, the free energy expression is relatively simple. In 
particular we find that the free energy of bending for a tubular membrane takes the approximate form, 

^F,,t-\kBTlp ds\{s)\ (11) 

where s is the the arc length and A(s) is the line curvature along the tube, while Ip is the persistence 
length of the tube: 

7rr(is:A + K + K||(l-V^)2 + K_L^2) 

= ksf • ^^^^ 

There are few experimental measurements of the persistence length of membrane tubes with protein 
lattices. For the F-BAR FBP17 producing tubes of radius r{FBP17) — 34nm the persistence length was 
measured to lp{FBP17) ~ 142/im [41] while for amphiphysin r(amph) ^ 7nm[51] and ^p(amph) = 9/im 
while for dynamin r(dynamin) ~ 20nm and Zp(dynamin) = 37/Ltm. A calculations of Ip from Eq.(12) 
solely based on k gives predictions which are an order of magnitude too small, which indicates that other 
elastic constants k\\,k± and gives the main contributions to Ip. 



E. Capping the tubes, Defects 



The formation of membrane tubes with helical coats seems to be generic for systems with membranes 
with membrane nematogens. Either the helical coat has to terminate resulting in an interfacial curve 
separating the coated and uncoated regions or the vesicle should sprout tubes and buds with the tips 
having a pair of point defects. The way this takes place in the tube end or at a domain boundary is 
mainly determined by the competition between interfacial tension, which in our model is largely regulated 
by the parameter J, and bending modulus. In Fig(7) is shown that when the interaction parameter J 
is increased the interfacial line shrinks, first transforming the vesicle from a disk to a structure with 
partially coated tubes and buds, but still no defects. Further increase in the line tension will result in 
tubes and buds that are fully coated but minimizing the length of the interfacial line between coated and 
uncoated regions. It does so by either moving the interfaces to the end of the tube forming a pair of 
point defects or deforming the membrane to form a narrow neck. Note that the line does not shrink to a 
single point defect of strength +1 but instead forms a pair of +1/2 defects bound to each other. This is 
the result of of the tt symmetry of the membrane nematogen and the strong coupling between membrane 
curvature and nematic orientation [52-54] . To our knowledge no details about the capping of the coated 
membrane tubes have been provided by experiments. 
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(a) J=0 (b) J=3 (c) J=5 

FIG. 7: A partly decorated membrane with k = 20, K|| = 5, ell = 3. Shown are: (a) a disc without a defect for 
4>A = 0.4, //y" = 0.5 and J = 0. (b) Tubes without defects at J = 3. ( the bottom panel shows an enlarged side 
view of the tip of a tube without defects ). (c) Tubes and buds with defects when J = 5. ( the bottom panel 
shows an enlarged top view of the tip of a bud with defects ) 



F. Curvature differences leads to segregation 




(a) (b) 

FIG. 8: (a) Tubular membrane of uniform cross section with fields A and B in the mixed phase for A) -^y* b — 
0.5, 0.5. (b) Activation of field B takes J/y* g from 0.5 to 1.25. The difference in the spontaneous curvatures leads 
to phase segregation and results in the constriction of tubes. Field concentrations are respectively i^a ~ 0.4 and 
0s = 0.6 while k = lOfcsT, ^ K.f ^ ^ksT and eaa = eiss = SksT. 

Another example of the curvature driven aggregation is demonstrated in Fig. 8. Shown in Fig. 8(a) 
is a tubular membrane of uniform cross section, fully decorated by two different types of membrane 
nematogens, labeUed A and B. The tube is stable, in the mixed state, when the directional spontaneous 
curvature of the in-plane fields, = b ^ ^-^^ ^'^'^ ^^"^ same. If a source of activation increases the 

spontaneous directional curvature of B to _ffO|| ^ — 1.25, analogous to activation of dynamin proteins 
by hydrolysis of GTP, the fields demix. The regions containing field B constrict the tube further. The 
equilibrium shape of the activated membrane is observed to have successive tubular regions of large 
and small curvatures ( see Fig. 8(b)), similar to the striated patterns of dynamin tubes obtained on 
treatment with GTP7S [55, 56]. For dynamin the molecular conformation and membrane tube diameter 
is GTP dependent [56, 57]. Furthermore it is observed that the tube constriction involves a tube twisting, 
suggesting a change in the helical angle [58]. In in vivo experiments, it has been demonstrated that 
structurally similar F-BAR proteins can co-localize into the same membrane tubes [14] while differing 
BAR proteins, like F-BAR and N-BAR, seggregate into separate membrane tubes with their respective 
characteristic r and Lp [41, 59]. Our analysis suggests that this recruitement of differing BAR-proteins 
into separate domains is possibly driven by their directional spontaneous curvatures. 
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G. Sheets versus tubes 

The effect of concentration, shown in Fig. (3), for surface coverages in the range 10% and 70% display 
a series of complex shape deformations connected to different aggregate structures. The regime, where 
inclusions stay separate, for the model parameters chosen here, appear at very low concentrations. The 
figure illustrates that for a system where the direct interactions parameter J between the membrane 
nematogens are weak the oligomers tend form larger rim-like formations, which stabilizes disc like struc- 
tures of vesicles. The rims form the edges of the discs. As the concentration is increased part of the edge 
turn tubular. 

So for a range of concentrations the disc and the tubes coexist. The tubules get more pronounced and 
the discs diminishes with the increase in concentration of membrane nematogens. Recent experiments 
on the formation of tubular (or smooth) ER suggest that some membrane curvature active proteins, 
reticulon protein and DP1[15], are highly enriched in the tubular ER [60] and the ER sheet edges[16]. 
Our results are thus in line with the idea that the concentration of these membrane nematogens are a 
major determinant for the amount of ER in sheets or tubules[16]. 

IV. CONCLUSION 

We have described the membrane curvature modifying properties of anisotropic protein inclusions, 
like the BAR proteins, in terms of an in-plane nematic field. We have shown that the flexibility of 
the membrane can promote aggregation and lateral domain formation of these membrane nematogens, 
even in the absence of self interactions. These domains can facilitate shape changes of the membrane. 
The equilibrium shapes obtained are strikingly similar to that seen in experiments involving curvature 
modifying proteins. Prominent structures seen are tubes and discs and coexistence of them. Depending 
on the preferred curvature of the nematogens, a protein lattice with helical nematic orientation around the 
tube is seen. The properties of this liquid crystalline structure was further analyzed from a continuum 
version of the model and the dependence of the tube radius and the orientation of the nematic with 
respect to the tube axis was calculated. We also estimate the thermally induced fluctuation in these 
quantities and show that they are comparable to what is seen in experiments. In addition we calculate 
the persistence length of the nematogen induced tubes and show that it is in the range of experimentally 
obtained values. This analysis provides the necessary basis to obtain estimates of model parameters from 
experiments on coated membrane tubes. At present the available experimental data are very limited. 

The present modeling provides additional support to the growing notion of the importance of local 
curvature modulating proteins in membrane shape generation in biological cells. Compared to previous 
modeling of the role of membrane proteins inducing directional membrane curvature we have taken into 
account that membranes are not fully decorated, the in-plane interactions between nematogens and the 
arbitrary membrane shapes with spherical topology. The current work focusses only on the membrane 
interacting part of the protein. Electric charges in BAR-protein are mostly localized to its membrane 
facing domain which in turn interact with anionic lipids and enables them to bind strongly to the 
membrane. One natural extension of this model is include the electrostatic interactions through cytosol 
between proteins moieties protruding out of the membrane. 

We emphasize that the main aim of this work is to show that anisotropic curvature induced by the 
inclusions can lead to aggregation and interesting shape changes. This is in contrast to the prevailing 
assumption that explicit protein interactions are essential for aggregation and formation of protein lattices. 
Though a quantitative comparison between the predictions of this model and experiments is not so easy, 
the model does demonstrate the possibility of generating many biologically relevant shapes of the vesicle 
by membrane mediated interactions alone. 
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